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Abstract 

The production of prompt charmed mesons D**, D+ and D*+, and their antiparticles, was measured 
with the ALICE detector in Pb-Pb collisions at the centre-of-mass energy per nucleon pair, 
of 2.76 TeV. The production yields for rapidity |y| < 0.5 are presented as a function of transverse 
momentum, pi, in the interval 1-36 GeV/c for the centrality class 0-10% and in the interval 1- 
16 GeV/c for the centrality class 30-50%. The nuclear modification factor /?aa was computed using a 
proton-proton reference at ^ = 2.76 TeV, based on measurements at -Js = 1 TeV and on theoretical 
calculations. A maximum suppression by a factor of 5-6 with respect to binary-scaled pp yields is 
observed for the most central collisions at px of about 10 GeV/c. A suppression by a factor of about 
2-3 persists at the highest pi covered by the measurements. At low pi (1-3 GeV/c), the /?aa has 
large uncertainties that span the range 0.35 (factor of about 3 suppression) to 1 (no suppression). In 
all px intervals, the /?aa is larger in the 30-50% centrality class compared to central collisions. The 
D-meson is also compared with that of charged pions and, at large p-x, charged hadrons, and 
with model calculations. 
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1 Introduction 


A state of strongly-interacting matter characterised by high energy density and temperature is predicted 
to be formed in ultra-relativistic collisions of heavy nuclei. According to calculations using Quantum 
Chromodynamics (QCD) on the lattice, these extreme conditions lead to the formation of a Quark-Gluon 
Plasma ((^P) state, in which quarks and gluons are deconfined, and chiral symmetry is partially restored 
(see e.g. 

Heavy quarks are produced in the hard scattering processes that occur in the early stage of the collision 
between partons of the incoming nuclei. Their production is characterised by a timescale A? < 1/ 

~ 0.1 fm/c for charm and ~ 0.01 fm/c for beauty quarks, that is shorter than the formation time of the 
QGP medium, about 0.3 fm/c at Large Hadron Collider (LHC) energies JH]. They can successively in¬ 
teract with the constituents of the medium and lose part of their energy, via inelastic processes (gluon 
radiation) |@, 0] or elastic scatterings (collisional processes) iM]. Energy loss can be studied using 
the nuclear modification factor /?aa> which compares the transverse-momentum (pj) differential produc¬ 
tion yields in nucleus-nucleus collisions (d^AA/dpr) with the cross section in proton-proton collisions 
(dapp/dpx) scaled by the average nuclear overlap function ((Taa)) 


^aa(pt) 


1 d^AA/dPT 
(Taa) dapp/dpT ' 


( 1 ) 


The average nuclear overlap function (Taa) over a centrality class is proportional to the number of 
binary nucleon-nucleon collisions per A-A collision in that class and it can be estimated via Glauber 
model calculations QQ. 

According to QCD calculations, quarks are expected to lose less energy than gluons because their cou¬ 
pling to the medium is smaller iaul]- In the energy regime of the LHC, light-flavour hadrons with px 
ranging from 5 to 20 GeV/c originate predominantly from gluons produced in hard scattering processes, 
while for larger px they originate mainly from light quarks (see e.g. 11311 1. Charmed mesons, instead, 
provide an experimental tag for a quark parent at all momenta. Therefore, the comparison of the heavy- 
flavour hadron /?aa with that of pions is expected to be sensitive to the colour-charge dependence of 
energy loss. However, other aspects than the energy loss, like the parton px spectrum and fragmentation 
into hadrons, influence the nuclear modification factor (see e.g. il3Lll4l1). The effect of the colour-charge 
dependence of the energy loss should be then studied via the comparison with model calculations, that 
include the description of the aforementioned aspects. 


In addition, several mass-dependent effects are predicted to influence the energy loss for quarks (see 111511 
for a recent review). The dead-cone effect should reduce small-angle gluon radiation for quarks that have 
moderate energy-over-mass values, i.e. for c and b quarks with momenta up to about 10 and 30 GeV/c, 
respectively 116142211 . Likewise, collisional energy loss is predicted to be reduced for heavier quarks, be¬ 
cause the spatial diffusion coefficient, which regulates the momentum transfers with the medium, scales 
with the inverse of the quark mass for a given quark momentum 12311 . In particular, the study of D mesons 
from low-px to high-px allows to study the variation of the energy loss for different charm quark veloc¬ 
ity: from a non-relativistic regime to an highly relativistic one. Low-momentum heavy quarks, including 
those shifted to low momentum by parton energy loss, could participate in the collective expansion of 
the system as a consequence of multiple interactions 1241. 12511. It was also suggested that low-momentum 
heavy quarks could hadronise not only via fragmentation in the vacuum, but also via the mechanism of 
recombination with other quarks from the medium 1251. 12611. 


The nuclear modification factor of heavy-flavour production was first studied at the Relativistic Heavy 
Ion Collider (RHIC). The PHENIX and STAR Collaborations reported measurements using heavy- 
flavour decay electrons and muons in Au-Au and Cu-Cu collisions at centre-of-mass energy per nu¬ 


cleon pair, = 200 GeV Il27l430ll . A suppression with respect to binary scaling was observed for px 
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larger than about 3 GeV/c, reaching a minimum Raa of about 0.2-0.3 in the interval 4 < pT < 6 GeV/c. 
The STAR Collaboration recently measured the Raa of D° mesons in Au-Au collisions for the interval 
0 < pt < 6 GeV/c 113IH . At pj of about 5-6 GeV/c the Raa value is similar to that observed for electrons 
from heavy-flavour decays and the Raa increases towards low pj, reaching a maximum value of about 
1.5 at 1-2 GeV/c. This feature is described by heavy-flavour transport calculations that include radial 
flow and a contribution due to recombination in the charm hadronisation process 113IH . 


A first measurement of the production of prompt D mesons at mid-rapidity in the pj interval 2-16 GeV /c 
was published, using the Pb-Pb data at = 2.76 TeV collected in 2010 during LHC Run 1 |32]. 
A minimum Raa of about 0.2-0.3 was measured at pj of about 10 GeV/c for the 20% most central 
collisions. The measurement of D-meson production in p-Pb collisions at = 5.02 TeV, showing an 
RpPb compatible with unity, has provided clear evidence that the suppression with respect to binary-scaled 
pp cross sections, observed in Pb-Pb collisions, cannot be attributed to cold nuclear matter effects for px 
larger than 2 GeV/c and is, thus, caused by final-state interactions in the hot and dense medium i33n . 


In Pb-Pb collisions, the nuclear modification factor at low px results from the interplay of different 
effects occurring in the initial and in the final state. The measured D-meson nuclear modification factor in 
p-Pb collisions, although consistent with unity, is also described within uncertainties by calculations that 
include substantial initial-state effects, such as parton shadowing or saturation [331], that could manifest 
as a reduction of the yields in Pb-Pb (and thus of the Raa) by up to 50% at low px. In addition, the 
measurement of a significant azimuthal anisotropy of D-meson production, with respect to the estimated 
direction of the reaction plane in non-central Pb-Pb collisions, indicates that chaim quarks participate 
in the collective expansion of the medium 041.13511. Therefore, radial flow could play a relevant role as 
well. In order to investigate these aspects, it is important to have a precise measurement of Raa down 
to low Px. In the high-px region, where parton energy loss is expected to be dominated by radiative 
processes, the extension of the D-meson Raa beyond 20 GeV/c would provide the first measurement of 
identified-hadron Raa at such high px. 


In this article we present the measurement of px-differential yields and nuclear modification factors of 
prompt D°, D+ and D*+ mesons (including their antiparticles), reconstructed via their hadronic decays 
in Pb-Pb collisions at = 2.76 TeV, using the data sample recorded in 2011. For central collisions, 
the integrated luminosity is larger by a factor of about 10 than that used for the previously published 
results [321] . This allows us to extend the measurement of Raa to lower and higher px (from 2-16 GeV/c 
to 1-36 GeV/c), to improve its precision, and to perform the study in a narrower class of the most central 
collisions (10% most central instead of 20% most central). 

The article is organised as follows: the experimental apparatus is described in Section|2l together with the 
data sample. In Section [3l the D-meson decay reconstruction and all corrections applied to the yields are 
presented, along with the procedure used to obtain the pp reference at = 2.76 TeV. In Section|4]the 
systematic uncertainties are discussed. The results for the 0-10% (central) and 30-50% (semi-peripheral) 
centrality classes are presented in Section [S] In the same Section results obtained in Pb-Pb collisions 
are compared with the nuclear modification factor measured in p-Pb collisions at y^^ = 5.02 TeV. 
A comparison with charged pions, charged particles (ch) and with theoretical model predictions is also 
reported. These comparisons are presented in terms of the ratio R^a/^aa well. Conclusions are 
drawn in Section [6] 


2 Experimental apparatus and data sample 


The ALICE experimental apparatus ll36ll is composed of various detectors for particle reconstruction 
and identification at mid-rapidity (Ifil < 0.9), a forward muon spectrometer (—4 < rj < —2.5) and a set 
of forward-backward detectors for triggering and event characterization. The detector performance for 
measurements in pp, p-Pb and Pb-Pb collisions from the LHC Run 1 is presented in 071] . 
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The main detector components used in this analysis are the VO detector, the Inner Tracking System (ITS), 
the Time Projection Chamber (TPC) and the Time Of Flight (TOE) detector, which are located inside a 
large solenoidal magnet providing a uniform magnetic held of 0.5 T parallel to the EHC beam direction 
(z axis in the AEICE reference system) and the Zero Degree Calorimeter (ZDC), located at ±114 m from 
the interaction point. 


Pb-Pb collision data were recorded with a minimum-bias interaction trigger based on information from 
the VO detector, which consists of two scintillator arrays covering the full azimuth in the pseudorapidity 
intervals —3.7<T]<—1.7 and 2.8 < T] < 5.1 li^ . The trigger logic required the coincidence of signals 
on both sides of the detector. An online selection based on the VO signal amplitudes was used to enhance 
the sample of central and mid-central collisions through two separate trigger classes. The scintillator 
arrays have an intrinsic time resolution better than 0.5 ns, and their timing information was used together 
with that from the ZDCs for offline rejection of events produced by the interaction of the beams with 
residual gas in the vacuum pipe. Only events with a reconstructed interaction point (primary vertex) 
within ±10 cm from the centre of the detector along the beam line were used in the analysis. 


Collisions were divided into centrality classes, determined from the sum of the VO signal amplitudes and 
defined in terms of percentiles of the total hadronic Pb-Pb cross section. In order to relate the centrality 
classes to the collision geometry, the distribution of the VO summed amplitudes was fitted with a function 
based on the Glauber model 111,12] combined with a two-component model for particle production S. 
The centrality classes used in the analysis are reported in Tab. [T] together with the average of the nuclear 
overlap function Taa, the number of events in each class (Vevents) and the integrated luminosity. 


The charged-particle tracks used to reconstruct the decay of D mesons were measured in the TPC and 
ITS. The tracking algorithm, based on a Kalman filter 14011 . starts from three-dimensional space points in 
the TPC, a large cylindrical drift detector with both total length and diameter of about 5 m, covering the 
pseudorapidity range |t]| < 0.9 with full azimuthal acceptance 14IH . Tracks are reconstructed in the TPC 
with up to 159 space points and with a measurement of the specific ionisation energy loss dE/dx with a 
resolution of about 6%. 


Hits in the ITS are associated to the prolongation of the TPC tracks, forming the global tracks. The ITS 
consists of six cylindrical layers of silicon detectors 14211 . The two innermost layers, placed at 3.9 and 
7.6 cm from the beam line, consist of Silicon Pixel Detectors (SPD). The third and fourth layers use 
Silicon Drift Detectors (SDD) and the two outermost layers contain double-sided Silicon Strip Detectors 
(SSD). The effective spatial resolutions, including the intrinsic detector resolutions and residual mis¬ 
alignments, are about 14, 40 and 25 /im, for SPD, SDD and SSD, respectively, along the most precise 
direction (r<p) 14211 . 


Global tracks are used to reconstruct the primary interaction vertex and the secondary vertices of D- 
meson decays. The transverse momentum resolution for global tracks ranges from about \% al pj = 
1 GeV/c to about 2% at 10 GeV/c, both in pp and Pb-Pb collisions. The spatial precision of global 
tracks is quantified by the resolution on the impact parameter do, which is the signed distance of closest 
approach between the track and the primary vertex in the vy-plane transverse to the beam direction. In 
Pb-Pb collisions, the do resolution is better than 65 pm for tracks with a transverse momentum larger 


Centrality class 

(Taa) (mb^^) 

^events 

^int (/^b ) 

0-10% 

23.44 ±0.76 

16.4 X 10^ 

21.3±0.7 

30-50% 

3.87±0.18 

9.0 X 10® 

5.8±0.2 


Table 1: Average of the nuclear overlap function, number of events and integrated luminosity for the two centrality 
classes used in the analysis. The uncertainty on the integrated luminosity stems from the uncertainty of the hadronic 
Pb-Pb cross section from the Glauber model ll39ll . 
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than 1 GeV/c and reaches 20 /xm for > 20 GeV/c 1371] . 

The TOE detector is an array of Multi-Gap Resistive Plate Chambers positioned at a distance of about 
370 cm from the beam line and covering the full azimuth over the pseudorapidity interval |t] | < 0.9. TOP 
particle identification is based on the difference between the particle arrival time at the TOP detector and 
a start time determined using the arrival time of all particles of the event with a minimization 1431] . 
The resolution (a) of the time-of-flight measurement is about 80 ps for pions at px = 1 GeV/c in the 
Pb-Pb collision centrality intervals used in this analysis. TOP provides charged-particle identification in 
the intermediate momentum range, with a 3a separation up to about 2.5 GeV/c for pions and kaons, and 
up to about 4 GeV/c for kaons and protons 1371] . 


3 Data analysis 


3.1 D-meson reconstruction 


DO, 

dO 


D+ and D*+ mesons, and their antiparticles, were reconstructed via their hadronic decay channels 
K^7r+ (weak decay with branching ratio, BR, of 3.88 ±0.05%), D+ —)■ K^;r+7r+ (weak decay, BR 


of 9.13 ±0.19%) and D*+ —D0;r+ (strong decay, BR of 67.7 ±0.05%) followed by D® —> K^;r+ i44l] . 
DO and D+ mesons have mean proper decay lengths (ct) of 123 and 312 pm, respectively 14411 . In the 
case of the D*+, the decay topology of the produced was exploited. The transverse momentum of the 
soft pions produced in the D*+ decays typically ranges from 0.1 to 1.5 GeV/c, depending on the D*+ 
Pt- 


and D+ candidates were formed using pairs and triplets of tracks with the correct charge-sign com¬ 
bination, requiring |t] | < 0.8, pj > 0.4 GeV/c, at least 70 associated space points (out of a maximum of 
159) and fit quality ;^^/ndf < 2 in the TPC, and at least two hits (out of six) in the ITS, out of which at 
least one in either of the two SPD layers. D*+ candidates were formed by combining candidates with 
tracks with px > 0.1 GeV/c and at least three hits in the ITS, out of which at least one in the SPD. 


The aforementioned track selection limits the D-meson acceptance in rapidity. The acceptance drops 
steeply to zero for |y| > 0.5 at low px and |y| > 0.8 for px > 5 GeV/c. A px-dependent fiducial ac¬ 
ceptance cut, IjdI < yfid(PT)> was therefore applied to the D-meson rapidity. The cut value, yfid(pT)> 
increases from 0.5 to 0.8 in the range 0 < px < 5 GeV/c according to a second-order polynomial func¬ 
tion and with a constant value of 0.8 for px > 5 GeV/c. 


The selection of the decay topology was based on the displacement of the decay tracks from the interac¬ 
tion vertex (via their impact parameter, do), the separation between the secondary and primary vertices 
(decay length, L) and the pointing angle of the reconstructed D-meson momentum to the primary vertex. 


This pointing condition was applied via a selection on the angle 


pointing 


between the direction of the 


reconstructed momentum of the candidate and the straight line connecting the primary and secondary 
vertices. The projections of the pointing angle and of the decay length onto the transverse plane (0j^jntj„g 
and L^) were also used. The selection requirements were tuned to provide a large statistical significance 
for the signal and to keep the selection efficiency as high as possible. The chosen selection values depend 
on the Px of the D meson and become ti ghte r from peripheral to central collisions. A detailed description 
of the selection criteria was reported in i32l .l35ll . 


In order to further reduce the combinatorial background, pions and kaons were identified using fhe TPC 
and TOP detectors. A 3a compatibility cut was applied to the difference between the measured and 
expected signals (for pions and kaons) for the TPC dE/dx and TOP time-of-flight. Tracks that are not 
matched with a hit in the TOP detector were identified using only the TPC information. Particle iden¬ 
tification (PID) was not applied to the pion track from the D*+ decay. This PID selection provides a 
reduction of the background by a factor of 2-3 at low px with respect to the case without applying the 
selection, while having an efficiency of about 95% for the signal. 
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Figure 1: (K, n) (top row) and (K, n, k) (central row) invariant-mass distributions for the centrality class 0-10%. 
Bottom row: Distribution of the mass difference AM ~ M{¥.nn) — M(K7r) for the centrality class 0-10%. The 
distributions are reported in three pj intervals for each meson (left, middle and right column). The fit functions 
used to describe the background (dash), the background without signal reflections (only for D**, long-dash) and the 
total distribution including the signal (solid) are shown. 


The raw D-meson yields were obtained from fits to the candidate invariant-mass distributions M(K;r) 
for D*’, M(K7r7r) for D+, and the mass difference AM = M(K7r7r) — M(K7r) for D*+. The D*’ and D+ 
candidate invariant-mass distributions were fitted with a function composed of a Gaussian for the signal 
and an exponential term to describe the background shape. In the 0-10% centrality class, the background 
in the M(K7r) distribution for the interval 1 < px < 2 GeV/c could not be accounted for by an exponential 
shape and was instead modelled with a fourth-order polynomial function. The AM distribution of 
candidates was fitted with a Gaussian function for the signal and a threshold function multiplied by an 
exponential for the background: _ 

In the case of D° mesons, an additional term was included in the fit function to account for the background 
from “reflections”, i.e. signal candidates that remain in the invariant-mass distribution when the (K, tt) 
mass hypotheses for the two decay tracks are swapped. A study of simulations showed that about 70% 
of these reflections are rejected by the PID selection, while the residual contribution results in a broad 
invariant-mass distribution, which can be described using a sum of two Gaussians. In order to account for 
the contribution of reflections in the data (2-5% at low pi, about 10% at high pT)^ a template consisting 
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M'(K7t) (GeV/c^) 



Figure 2: (K, n) invariant-mass distribution for the interval I < pj <2 GeV /c for the 0-10% (left) and 30-50% 
(right) centrality classes, obtained after the subtraction of the background estimated by a fourth-order polynomial 
function for the most central collisions and an exponential for the 30-50% centrality class. The contribution of 
reflections is also included in the fit. The fit function used to describe the signal (solid line) is shown. 

of two Gaussians was included in the fit. The centroids and widths, as well as the ratios of the integrals of 
these Gaussians to the signal integral, were fixed to the values obtained in the simulations (see also Is^l ). 

In the centrality class 0-10%, the signal extraction was performed in the interval 1 < pj < 24 GeV/c 
for D° mesons, divided in 9 pt bins, and in the interval 3 < pj < 36 GeV/c for D+ and D*+ mesons, 
divided in 8 pj bins. In the centrality class 30-50%, the signal extraction was possible in the interval 

1 < Pt < 16 GeV/c for mesons and in the interval 2 < px < 16 GeV/c for D+ and D*+ mesons. 
Beyond these intervals, the signal extraction was prevented by the low signal-over-background ratio at 
low Px, and by the low signal yield at high px- Figure[T]shows the D° and D+ invariant-mass distributions 
and D*+ mass difference distributions in three px intervals for the centrality class 0-10%. In the interval 
16 < Px < 24 GeV/c the fit range for the D° case is asymmetric. The range was limited to values larger 
than 1.68 GeV/c^ because the invariant-mass distribution of (K, tt) pairs from D mesons decaying in 
three or more prongs produces a wide structure below about 1.72 GeV/c^, which cannot be accounted 
for by the background terms of the fit function. 

Figure |2] shows the background-subtracted D° invariant-mass distribution for the interval 1 < px < 

2 GeV/c for the 0-10% (left panel) and 30-50% (right panel) centrality classes. 

For all three D-meson species, the position of the signal peak was found to be compatible with the world 
average value and its px-dependent width with the values observed in the simulation. The statistical 
significance of the observed signals S/VS + B varies from 3 to 18, while the signal-over-background 
ratio S/B ranges from 0.01 to 1.8, depending on the meson species, px interval and centrality class. 

3.2 dA/dpx spectra corrections 

The D-meson raw yields were corrected in order to obtain the px-differential yields of prompt D mesons 


dA° 


dpx 


|y|<0.5 


/prompt(PT) • 


bl<yfid 


Apx ■ 0(v’ ' X S) prompt (Pt) ‘ BR • Aevents 


( 2 ) 
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where prompt refers to mesons not coming from weak decays of B hadrons. The raw yields were 

divided by a factor of two to obtain the charge-averaged (particle and antiparticle) yields. To correct for 
the contribution of B-meson decay feed-down, the raw yields were multiplied by the fraction of promptly 
produced D mesons, /prompt (discussed in details later in this section). Furthermore, they were divided 
by the product of prompt D-meson acceptance and efficiency (Ace x e)prompt! by the decay channel 
branching ratio (BR), by the transverse momentum interval width (Apx) and by the number of events 
(Aevents)- The factor cty. = yfid/0.5 normalises the corrected yields measured in |y| < y^d to one unit of 
rapidity |y| < 0.5, assuming a uniform rapidity distribution for D mesons in the measured range. This 
assumption was validated to the 1% level with simulations 


The correction for acceptance and efficiency (Ace x e)prompt was defermined using Monfe Carlo sim¬ 
ulations wifh a defailed description of the detector and its response, based on the GEANT3 transport 
package 14711 . The simulation was tuned to reproduce the (time-dependent) position and width of the 
interaction vertex distribution, as well as the number of active read-out channels and the accuracy of the 
detector calibration. The underlying Pb-Pb events at = 2.76 TeV were simulated using the HIJING 
vl.383 generat or and D-meson signals were added with the PYTHIA v6.421 generator with 
Perugia-0 tune 15011 . Each simulated PYTHIA pp event contained a cc or bb pair, and D mesons were 
forced to decay in the hadronic channels of interest for the analysis. Out of all the particles produced in 
these PYTHIA pp events, only the heavy-flavour decay producfs were kepi and Iransporled Ihrough Ihe 
defector simulafion fogefher with the particles produced according to HIJING. In order to minimise the 
bias on the detector occupancy, the number of D mesons injected into each HIJING event was adjusted 
according to the Pb-Pb collision centrality. In the most central event class, the px distribution of D 
mesons was weighted in order to match the shape measured for the D° meson. In the semi-peripheral 
centrality class, the D-meson px distribution was weighted so as to match the shape given by fixed- 
order-next-to-leading-log perturbative OCD calculations (FONEE) 115iL 15211 multiplied by the Raa{pt) 
computed using the BAMPS model 153145511 . 


The efficiencies were evaluated from simulaled evenfs lhal have Ihe same average charged-parlicle mulli- 
plicify, corresponding fo Ihe same defector occupancy, as observed in dafa in the centrality classes 0-10% 
and 30-50%. Figure|3]shows the D°, D+ and D*+ acceptance-times-efficiency (Ace x e) for primary and 
feed-down D mesons with rapidity |y| < yfid(pT) in the centrality class 0-10%. The efficiencies range 
from aboul 0.1% al low px to 10-30% at high px, because of the momentum dependence of the D-meson 
decay length and of the topological selections applied in the different momentum intervals. Also shown 
in the figure are fhe (Ace x e) values for the case where no PID is applied. The relative difference with 
respect to the (Ace x e) obtained using the PID selection is about 5%, illustrating the high efficiency of 
fhe PID criteria. The (Ace x e) for D mesons from B-meson decays is larger than for prompt D mesons 
by a factor of about 1.5, because the decay vertices of the feed-down D mesons are more separated from 
the primary vertex and are, therefore, more efficiently selected by the analysis cuts. 


The /prompt factor was obtained, following the procedure introduced in Il32ll . as 


/prompt — 1 


N, 


D+D feed-down 


D+D 


= l-R 


N 


feed-down 

AA 


(Taa)' 


FONLL.EvtGen 


da Y 

tlPT / feed-down, h|<0.5 


Apx * OCy * (Ace X £)feed-down * BR * Aei 

i/uD+D 
2 ^’'raw 


(3) 


In this expression, the symbols denoting the px dependence have been omitted for brevity, Ar'^/^ is the 
measured raw yields and is the estimated raw yields of D mesons from B-meson decays. In 

detail, the B-meson production cross section in pp collisions at = 2.76 TeV, estimated with FONEE 
calculations llsdi . was folded with the B —)• D -|-A decay kinematics using the EvtGen package |57] and 
multiplied by (Taa) in each centrality class, by the (Ace x e) for feed-down D mesons, and by the other 
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Figure 3: Product of acceptance and efficiency for D mesons in Pb-Pb collisions for the 0-10% centrality class. 
The rapidity interval is |y| < ygd (see Section 3.1). The values for prompt (solid lines) and feed-down (dotted lines) 
D mesons are shown. Also displayed, for comparison, are the values for prompt D mesons without PID selection 
(dashed lines). 


factors introduced in Eq. In addition, the nuclear modification factor of D mesons from B-meson 
decays was accounted for. The comparison of the 7 ?aa of prompt D mesons (7?^™™^^) 11581] with that 
of J/t/r from B-meson decays measured in the CMS experiment indicates that chaimed hadrons 


are more suppressed than beauty hadrons. The value _ 2 . was used to compute the 

?feed-down Z/^P^^P^ < 3 was Considered for the evaluation 


correction, and the variation over the range 1 < 7?^^ 


of the systematic uncertainties, in order to take into account possible centrality and px dependences. 


Assuming R 


feed-down 

AA 


= 2R 


,prompt 
AA 


, the resulting /prompt ranges from about 0.65 to 0.85, depending on the 


D-meson species and on the px interval. 

3.3 Proton-proton reference for 7?aa 

The px-differential cross section of prompt D mesons with |y| < 0.5 in pp collisions at ^/s = 2.76 TeV, 
used as reference for the nuclear modification factor, was obtained as follows: 


in the interval 2 < px < 16 (24) GeV/c for D° (D+ and D*+), the measurement at y/s = 7 TeV ll45h 
scaled to y/s = 2.76 TeV with EONEE calculations lli^ was used; 

in the interval 1 < px < 2 GeV/c for D*’, an average of the aforementioned yfs = 7 TeV scaled 
measurement and of the measurement at yfs = 2.76 TeV 0] was used; 

in the interval 16(24) < px < 24(36) GeV/c for D° (D+ and D*+), where their cross sections 
were not measured in pp collisions, the EONEE calculation at yfs = 2.76 TeV 11561] was used as a 
reference, after scaling it to match the central value of the data at lower px- 


The px-dependent scaling factor from yfs = 7 TeV to yfs = 2.76 TeV was determined with EONEE 
calculations and its uncertainties were determined by varying the parameters (charm-quark mass, fac¬ 
torisation and renormalisation scales) as described in 116IH . The uncertainties on the scaling factor range 
from for 1 < Px < 2 GeV/c to about ±5% for px > 10 GeV/c. The result of the scaling of the 
yfs = 7 TeV Px-differential cross sections to yfs = 2.76 TeV was validated with measurements from a 
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smaller data sample in pp collisions at y/s = 2.76 TeV 1^ . These measurements cover a reduced pj 
interval with a statistical uncertainty of 20-25% and therefore they were not used as a pp reference for 
Pj >2 GeV / c. 

For the lowest pj interval for the D*’ meson, the two references (obtained from the measurement in pp col¬ 
lisions at y/s = 2.76 TeV and from the y/s = 1 TeV scaled measurement) have comparable uncertainties. 
Therefore, in this interval, the two values were averaged using the inverse of the squared relative uncer¬ 
tainties as weights. In particular, the statistical uncertainties and the uncorrelated part of the systematic 
uncertainties, i.e. the systematic uncertainty from data analysis (yield extraction, efficiency corrections) 
and the scaling uncertainty, were used in the weight. The uncertainties on the feed-down subtraction 
were considered as fully correlated among the two measurements, and were propagated linearly. 


The cross section measurements for D mesons in pp collisions at yfs = 7 TeV are limited to px < 
16 GeV/c for and to pj < 24 GeV/c for and D*^. Beyond these lim its t he pp reference was ob¬ 
tained using the cross section from the FONLL calculation at y/s = 2.76 TeV 1^ . Since the central value 
of the FONFF calculation underestimates the measurement for pj > 5 GeV/c at both y/s = 2.76 TeV 


and y/s = 1 TeV i45l . 16011 . the FONFF cross section was multiplied by a scaling factor (fc) 


da 

dpT 


= K- 



(4) 


V7=2.76TeV,h|<0.5 


The factor fc was determined by fitting with a constant the data-to-theory ratio at y/s = 7 TeV in the 
interval 5 < /jt < 16 GeV/c. Since the measurements at y/s = 2.76 TeV are less precise, they do not 
constrain further the scaling factor. Depending on the D-meson species, the factor K ranges from 1.4 to 
1.5, for the central values of the FONFF calculation parameters 15611 . The statistical uncertainty of the 
extrapolated cross section was determined by propagating the statistical uncertainties of the measurement 
in the determination of K and it amounts to about 5%. The systematic uncertainties were evaluated under 
the conservative assumption that the systematic uncertainties of the measurement are fully correlated 
over pt, i.e. by repeating the calculation of K after shifting all data points consistently within their 
systematic uncertainties. In addition, the calculation in Eq. Q was performed considering the FONFF 
cross section s obtained from combinations of the renormalisation and factorisation scales with values 
(0.5,1,2) • + Px c 156], as well as the upper and lower limits of their envelopeQ. This resulted in a 


total systematic uncertainty on the px-extrapolated cross section of about -35%- 


4 Systematic uncertainties 

4.1 Systematic uncertainties on the D-meson pj spectra 

The systematic uncertainties were estimated as a function of transverse momentum for the two centrality 
classes. Table|2]lists the uncertainties for three pj intervals for each meson species. 

The systematic uncertainty on the raw yield extraction was evaluated by repeating the fit of the invariant- 
mass distributions while varying the fit range; by fixing fhe mean and sigma of fhe Gaussian ferm fo 
fhe world-average value and fhe expecfafions from Monfe Carlo simulations, respectively; and by using 
differenf fif functions for fhe background. Specifically, firsf- and second-order polynomials were used 
for D° and D+, and a power law mulfiplied by an exponenfial or a fhreshold funcfion for D*+. A mefhod 
based on bin counfing of the signal after background subtraction was also used. This method does not 
assume any particular shape for the invariant-mass distribution of the signal. The estimated uncertainties 
depend on the centrality class and on the pj interval, ranging from 5% to 15% for D*^, 8% to 10% for D+ 
and 5% to 10% for D*+, typically with larger values in the lowest and highest pj intervals. 

'where me and px.c are respectively the mass and the transverse momentum of the charm quark considered in the calcula¬ 
tions. 
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Particle 

D° 

D+ 

D*+ 

0-10% centrality class 

Px interval (GeV/c) 

1-2 

6-8 

16-24 

3-4 

6-8 

24-36 

3-4 

6-8 

24-36 

Yield extraction 

15% 

5% 

15% 

10% 

8% 

8% 

12% 

5% 

10% 

Tracking efficiency 

10% 

10% 

10% 

15% 

15% 

15% 

15% 

15% 

15% 

Selection cuts 

15% 

5% 

5% 

10% 

10% 

10% 

10% 

10% 

10% 

PID efficiency 

5% 

5% 

5% 

5% 

5% 

5% 

5% 

5% 

5% 

MC Px shape 

15% 

1% 

1% 

6% 

1% 

1% 

4% 

1% 

1% 

FONLL feed-down corr. 

^ 45 % 

-13 

t\i% 

ii2% 

-11 

+ 14 % 

+ 3% 
-12 


1 8% 

^fed-down 

+ 5% 

+ 11% 
-10 

-> 

+ ^% 

+ ^% 




- 6 % 

BR 


1.3% 



2.1% 



1.5% 


Centrality class definition 


< 1% 



< 1% 



< 1% 


30-50% centrality class 

Px interval (GeV / c) 

1-2 

6-8 

12-16 

2-3 

6-8 

12-16 

2-3 

6-8 

12-16 

Yield extraction 

10% 

8% 

8% 

10% 

10% 

12% 

12% 

8% 

5% 

Tracking efficiency 

10% 

10% 

10% 

15% 

15% 

15% 

15% 

15% 

15% 

Selection cuts 

10% 

10% 

15% 

10% 

10% 

15% 

15% 

10% 

5% 

PID efficiency 

5% 

5% 

5% 

5% 

5% 

5% 

5% 

5% 

5% 

MC Px shape 

5% 

1% 

3% 

10% 

2% 

2% 

10% 

1% 

1% 

FONFF feed-down corr. 

^ 45 % 

^ 12 % 

+ 

-11 

^ 2 > 

^ 12 % 

+“% 

-13 

+ 3 % 
-19 

^ 8% 

+ > 

^fe^d-down/^mt 

+ 6% 


-11 

+ 1 % 


+ 16 % 

-12 


+ ’% 

1^% 

BR 


1.3% 



2.1% 



1.5% 


Centrality class definition 


2% 



2% 



2% 



Table 2: Relative systematic uncertainties on the prompt D-meson production yields in Pb-Pb collisions for three 
selected pj intervals, in the two centrality classes. 


For D° mesons, the systematic uncertainty due to signal reflections in the invariant-mass distribution was 
estimated by changing by zb50% the ratio of the integral of the reflections over the integral of the signal 
(obtained from the simulation) used in the invariant-mass fit with the reflections template. In addition, 
the shape of the template was varied using a polynomial parameterisation (of third or sixth order) of the 
simulated distribution, instead of a double-Gaussian parameterisation. A test was carried out using, in the 
fit, a template histogram of the reflections obtained directly from the simulation, rather than a functional 
form. The variation observed in the raw yields, ranging from 3% to 7% from low to high pi, was added 
in quadrature as an independent contribution to the yield extraction systematic uncertainty. 


The systematic uncertainty on the tracking efficiency was estimated by comparing the probability to 
match the TPC tracks to the ITS hits in data and simulation, and by varying the track quality selection 
criteria (for example, the minimum number of associated hits in the TPC and in the ITS, and the maxi¬ 
mum ;^^/ndf of the momentum fit). The efficiency of the track matching and the association of hits in the 
silicon pixel layers was found to be well reproduced by the simulation with maximal deviations on the 
level of 5% in the pj range relevant for this analysis (0.5-25 GeV/c) 1371]. The effect of mis-associating 
ITS hits to tracks was studied using simulations. The mis-association probability is about 5%, for central 
collisions, in the transverse momentum interval 1 < px < 3 GeV /c and drops rapidly to zero at larger px. 


11 




















Transverse momentum dependence of D-meson production in Pb-Pb collisions ALICE Collaboration 


It was verified that the signal selection efficiencies are the same for D mesons with and without wrong 
hit associations. The total systematic uncertainty on the track reconstruction procedure amounts to 5% 
for single tracks, which results in a 10% uncertainty for D*’ mesons (two-tracks decay) and 15% for D+ 
and D*+ mesons (three-tracks decay). 

The uncertainty on the D-meson selection efficiency reflects a possible non-exact description of the D- 
meson kinematic properties and of the detector resolutions and alignments in the simulation. This effect 
was estimated by repeating the analysis with different values of the selection cuts, significantly modifying 
the efficiencies, raw yield and background values. As expected, larger deviations in the corrected yields 
were observed at low pi, where the efficiencies are low and vary steeply with ptj because of the tighter 
selections. Due to this, the systematic uncertainties are slightly larger in these px intervals. The assigned 
systematic uncertainty varies from 5% to 15% for D*’, equals 10% for D+, and varies from 10% to 15% 
for D*+. 


A 5% systematic uncertainty related to the PID selection was evaluated by comparing the ratio of the 
corrected yields extracted with and without particle identification. 

The uncertainty on the efficiencies arising from the difference between the real and simulated D-meson 
transverse momentum distributions depends on the width of the px intervals and on the variation of the 
efficiencies within them. This uncertainty also includes the effect of the px dependence of the nuclear 
modification factor. As explained in Sectionfor the centrality class 0-10%, the D-meson transverse 
momentum distribution from the PYTHIA simulation was re-weighted in order to reproduce the D° 
spectrum shape observed in data, while for the 30-50% centrality class, the weights were defined in 
order fo match the px distributions from FONLL calculations multiplied by the /?aa from the DAMPS 
model. A systematic uncertainty was estimated by using two alternative D-meson px distributions in 
both centrality classes: i) FONLL px distributions, ii) FONLL px distributions multiplied by 7 ?aa from 
the DAMPS model. In addition, for the most central events, a different parameterisation of the measured 
Px spectrum was used. The resulting uncertainties decrease with increasing px, varying from 5-6% to 
1% in the interval 2 < px < 36 GeV/c. For D*’ mesons, efficiencies increase by more than a factor 
five wifhin the interval 1 < px < 2 GeV/c in the most central collisions. As a consequence, a larger 
uncertainty of 15% resulted from a detailed study of the stability of the corrected yields when changing 
the Px spectrum in the simulation. 


The systematic uncertainty on the subtraction of feed-down from D decays (i.e. the calculation of the 
/prompt fraction) was estimated i) by varying the px-differential feed-down D-meson cross section from 
the FONLL calculation within the theoretical uncertainties, ii) by varying the hypothesis on the ratio of 
the prompt and feed-down D-meson /?aa in the range 1 < < 3, and iii) by applying 

an alternative method to compute /prompt- This second method is based on the ratio of charm and beauty 
FONLL cross sections, instead of the absolute beauty cross section. The procedure is the same used 
for previous measurements of D-meson production with ALICE ||^, 35, i3l- The resulting uncertainty 
ranges between ^ 45 % at low px and /g% at high px for the 0 - 10 % centrality class, and between ^ 45 % 
at low Px and /g% at high px for the 30-50% centrality class. The uncertainty from the variation of 
the feed-down D-meson /?aa hypothesis ranges from 6 to 16%, as shown in Fig. |4j where the relative 
variation of the prompt D° yield is shown as a function of the hypothesis on for four 

Px intervals. 


The uncertainties on the branching ratios were also considered ll44ll as well as the contribution due to 
the 1.1% uncertainty on the fraction of the hadronic cross section used in the Glauber fit to determine 
the centrality classes 11391] . The latter was estimated from the variation of the D-meson dA/dpx when 
the limits of the centrality classes are shifted by ±1.1% (e.g. shifted from 30-50% to 30.3-50.6% and 
29.7-49.5%) 1321]. The resulting uncertainty, common to all px intervals, is smaller than 1% for the 
0-10% centrality class and about 2% for the 30-50% centrality class. 
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Hypothesis on feed-down)/(R^_i^ prompt) 


Figure 4: Relative variation of the prompt D** yields as a function of the hypothesis on down^^piompt 
B-meson feed-down subtraction. 


4.2 Systematic uncertainties on /?aa 

The systematic uncertainties on the Raa measurement include those on the D-meson corrected yields, 
those on the proton-proton cross section reference, and the uncertainties on the average nuclear overlap 
function. 


The systematic uncertainties on the D-meson corrected yields are obtained considering as uncorrelated 
the different contributions described in the previous section. 


The uncertainty on the pp reference used for the calculation of Raa has two contributions. The first 
is the systematic uncertainty on the measured px-differential D-meson cross section at ^/s = 1 TeV. 
This uncertainty is about 25% at the lowest pi and 17% at the highest px for D*’ mesons, excluding 
the uncertainty for feed-down corrections, and few percent larger for D+ and D*+ mesons idSH . The 
systematic uncertainty on the feed-down subtraction deriving from the variation of the parameters of the 
FONLL calculation and from the use of the alternative method to compute /prompt was considered to be 
correlated in the Pb-Pb and pp measurements. These variations were carried out simultaneously for the 
numerator and denominator of Rak, so only the residual effect was attributed as a systematic uncertainty. 
Therefore, the variation of the value of 
the feed-down uncertainty on /?aa- 


.prompt 2 3 jg {j^g main contribution to 


The second contribution to the pp reference uncertainty is the scaling to = 2.76 TeV. It ranges from 
/jQ% in the interval 2 < px < 3 GeV jc to about 5% for px > 10 GeV/c |^]. Note that the upper/lower 
uncertainties are reversed when considering /?aa> where the pp reference is in the denominator. In the 
interval 1-2 GeV/c, this scaling uncertainty is much larger (/jj%), but its impact on the pp reference 
was reduced by about a factor of two by using a weighted average of the cross section scaled from 7 TeV 
and the measured cross section at 2.76 TeV (see Section lT^ . 


The extrapolation of the pp reference to the intervals 16 < px < 24 GeV/c for D*^ mesons and 24 < px < 
36 GeV/c for D+ and D*+ mesons resulted in a total systematic uncertainty of about ^ 59 %, as described 
in Section [TH 
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Particle 

D® 

D+ 

D*+ 

0 -10% centrality class 

Px interval (GeV/c) 

1-2 

6-8 

16-24 

3-4 

6-8 

24-36 

3^ 

6-8 

24-36 

dApb-Pb/dpx (excl. feed-down) 

28% 

14% 

22 % 

22 % 

20 % 

22 % 

24% 

20 % 

21 % 

dApp/dpx (excl. feed-down) 

21 %* 

16% 

17% 

20 % 

19% 

20 % 

17% 

17% 

18% 

— scaling of the pp ref. 

13^%* 

^10% 

- 

+ 

-19 

+ 6% 
-10 

- 

+ 9% 
-20 '® 

-> 

- 

High-px extrapolation 

- 

- 

+34% 

-51 

- 

- 

+37% 

-56'® 

- 

- 

^53% 

FONLL feed-down corr. 

1 1% 

^ ^% 

+ 10% 
-16'® 


1 1% 

+ 

-14 '® 


t 1% 

^ 8% 


1^9% 

-11 

-13 

t ?% 

1 ^2% 

+ 16% 
-12 '® 


+ ^% 

+ ^% 

Normalisation 


4.8% 



4.8% 



4.8% 


30-50% centrality class 

Px interval (GeV/c) 

1-2 

6-8 

12-16 

2-3 

6-8 

12-16 

2-3 

6-8 

12-16 

dApb-pb/dpx (excl. feed-down) 

20 % 

20 % 

22 % 

25% 

21 % 

22 % 

29% 

19% 

18% 

dApp/dpx (excl. feed-down) 

21 %* 

16% 

17% 

20 % 

19% 

20 % 

17% 

17% 

18% 

— scaling of the pp ref. 

13^%* 

^10% 


-19 m 


^6% 

+ 9% 
-20 '® 

^1^® 

1^% 

FONLL feed-down corr. 


^3% 


t 2% 

+ 1% 



+ 3% 

^ 3% 

^feed-down ^^prompt 

1^9% 

-11 

+ 15% 
-11 

1 ?% 

+13% 
-10 '® 

-13 '® 

-1% 

1 + 

oo o 


Normalisation 


6 .2% 



6 .2% 



6 .2% 



Table 3: Relative systematic uncertainties on the prompt D-meson Raa for three px intervals, in the two centrality 
classes. Uncertainties marked with a * were obtained as the average of the measurement at i/s = 2.76 TeV and the 
measurement at = 1 TeV, scaled using FONLL if^ . as described in Section iTJ] 


The uncertainties on /?aa are listed in Tab. |3] The uncertainties on the normalisation are the quadratic 
sum of the pp normalisation uncertainty (3.5%) and the uncertainty on {Taa), which is 3.2% and 4.7% 
in the 0-10% and 30-50% centrality classes, respectively. 

All the uncertainties described in this Section that result from detector effects are considered to be largely 
correlated over transverse momentum, with the exception of the yield extraction uncertainty that depends 
on the S/B in each pj interval. The uncertainties related to the feed-down assumptions and to the ^/s- 
scaled pp reference are fully correlated over pj, with the exception of that for the hypothesis on the ratio 
of the prompt and feed-down D-meson Raa that might not be constant as a function of pj. 

5 Results and discussion 
5.1 D-meson pj spectra and Raa 

The transverse momentum distributions dN/dpj of prompt D^, D+ and D*+ mesons are shown in 
Figs. |5(a)| |5(b)| and |5(c)| for the 10% most central Pb-Pb collisions. The results are presented in the 
interval 1 < pj < 24 GeV/c for the D° mesons and 3 < pj < 36 GeV/c for D+ and D*+ mesons. They 
are compared to the corresponding pp cross section reference multiplied by (Taa)- The vertical bars rep¬ 
resent the statistical uncertainties, the empty boxes the systematic uncertainties from the data analysis, 
and the shaded boxes the systematic uncertainty due to the subtraction of the feed-down from B-hadron 
decays. Uncertainties on the pp cross section normalisation and on the branching ratios are quoted sep¬ 
arately. A clear suppression of the D-meson yields is observed at intermediate (3 < px < 8 GeV/c) and 
high transverse momenta (px > 8 GeV/c) in central Pb-Pb collisions as compared to the binary-scaled 
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Figure 5: Transverse momentum distributions dN/dpj of prompt D** (a), D+ (b) and D*+ (c) mesons in the 0-10% 
centrality class in Pb-Pb collisions at ^Ann = 2.76 TeV. The pp reference distributions ( 7 aa) dcr/dpT are shown 
as well. Statistical uncertainties (bars) and systematic uncertainties from data analysis (empty boxes) and from 
feed-down subtraction (shaded boxes) are shown. Horizontal bars represent bin widths, symbols are placed at the 
centre of the bin. The dN/dpj distributions of the three D-meson species in the 10% most central Pb-Pb collisions 
are compared to each other in panel (d), where the D*+ production yields are scaled by a factor of five for visibility. 


pp reference. In Fig. |5(d)| the transverse momentum distributions of prompt D°, D+ and D*+ mesons 
in the 10% most central collisions are compared to each other. The dN/dpj values of D*+ mesons are 
scaled by a factor of five for visibility. 


The D-meson dN/dpj distributions measured in the 30-50% centrality class are shown in Fig.|6] Also 
for this centrality class, a clear suppression of the D-meson yields as compared to the expectation based 
on binary scaling of the pp yields is observed for pt > 3 GeV/c. In Fig. |6(d)[ the dN /dpj of prompt D°, 
D+ and D*+ (the latter scaled by a factor of five) are compared to each other. 
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Figure 6: Transverse momentum distributions dN/dpj of prompt D** (a), D+ (b) and D*+ (c) mesons in the 30- 
50% centrality class in Pb-Pb collisions at = 2.76 TeV. The pp reference distributions (7aa) da/dpj are 
shown as well. Statistical uncertainties (bars) and systematic uncertainties from data analysis (empty boxes) and 
from feed-down subtraction (shaded boxes) are shown. Horizontal bars represent bin widths, symbols are placed 
at the centre of the bin. The dN/dpi distributions of the three D-meson species in Pb-Pb collisions in the 30-50% 
centrality class are compared to each other in panel (d), where the D*+ production yields are scaled by a factor of 
five for visibility. 
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(GeV/c) 



p^ (GeV/c) 


Figure 7: Ratio of prompt D-meson yields (D+ZD*^ and D*+/D‘') as a function of pj in the 10% most central 
Pb-Pb collisions at = 2.76 TeV compared to the results in pp collisions at yG = 7 TeV. Statistical (bars) and 
systematic (boxes) uncertainties are shown. 



Figure 8: Raa of prompt D**, D+, and D*+ mesons for the 0-10% (left) and 30-50% (right) centrality classes. 
Statistical (bars), systematic (empty boxes), and normalisation (shaded box) uncertainties are shown. Horizontal 
bars represent bin widths. D** symbols are placed at the centre of the bin. D+, and D*+ are shifted for visibility. 


Figure |7] shows the pT-dependent ratios of D+/D° and D*+/D® for central Pb-Pb collisions. They are 
found to be compatible within uncertainties with those measured in pp collisions at ^/s = 7 TeV 16211 . 
Similar results were also found for the 30-50% centrality class. Therefore, no modification of the relative 
abundances of these three D-meson species is observed within the current uncertainties in central and 
semi-central Pb-Pb collisions relative to the pp ones at LHC energies. 


The Raa of prompt D*’, D+ and D*+ mesons is shown in Fig. [8] for the 0-10% (left panel) and 30- 
50% (right panel) centrality classes. The statistical uncertainties, represented by the vertical error bars, 
range from 10% in the intermediate pj range up to about 25-30% in the lowest and highest pj inter¬ 
vals, for the 10% most central collisions. The statistical uncertainty on the reference measurement at 
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y/s = 7 TeV dominates this uncertainty in the interval 2 < pj < 16 GeV/c. For the 30-50% centrality 
class, the statistical uncertainties at low and intermediate pj are similar in magnitude to those of central 
collisions and are about 20% in the interval 12 < pj < 16 GeV/c. The total px-dependent systematic 
uncertainties, described in the previous Section, are shown as empty boxes. The normalisation uncer¬ 
tainty is represented by a filled box at /?aa = 1 ■ The nuclear modification factors of the three D-meson 
species are compatible within statistical uncertainties for both centrality classes. For the 10% most cen¬ 
tral collisions, the measured Raa shows a suppression that is maximal at around pj = 10 GeV/c, where 
a reduction of the yields by a factor of 5-6 with respect to the binary-scaled pp reference is observed. 
The suppression decreases with decreasing px for px < 10 GeV/c, and it is of the order of a factor of 3 
in the interval 3 < px < 4 GeV/c, while the Raa ranges from about 0.35 to 1 in the first two px intervals. 
For Px > 10 GeV/c, the suppression appears to decrease with increasing px, but the large statistical 
uncertainties do not allow us to determine the trend of the Raa- A suppression (Raa < 0.5) is still ob¬ 
served for D mesons with px > 25 GeV/c. For the 30-50% centrality class, the suppression amounts 
to about a factor of 3 at px = 10 GeV/c, which indicates that the suppression of the high-px D-meson 
yields is smaller than in the 0-10% centrality class. As for the central collisions, the suppression reduces 
at lower momenta, with Raa increasing with decreasing px up to a value of about 0.6 in the interval 
3 < Px < 4 GeV/c. For lower px the suppression is further reduced and Raa is compatible with unity. 


The average nuclear modification factor of D°, D+ and D*+ mesons was computed using the inverse 
of the squared relative statistical uncertainties as weights. The systematic uncertainties were propagated 
through the averaging procedure, considering the contributions from the tracking efficiency, the B-meson 
feed-down subtraction and the FONFF-based ^^-scaling of the pp cross section from y/s = 7 TeV to 
y/s = 2.76 TeV as fully correlated uncertainties among the three D-meson species. The average D- 
meson Raa for the two centrality classes is shown in the left panel of Fig. |9l A larger suppression, by 
about a factor of two, is observed in the 10% most central collisions compared to the 30-50% centrality 
class for px > 5 GeV/c. The stronger suppression observed in central collisions can be understood as 
resulting from to the increasing medium density, size and lifetime from peripheral to central collisions. 
The Raa values measured for the 0-10% centrality class are slightly lower, although compatible within 
uncertainties, than those reported in Ref. 0211 for the 20% most central collisions, measured with the 
2010 data sample. As a consistency check, the analysis on the 2011 data sample was also performed 
in the 0 - 20 % centrality class and the resulting Raa value was found to be compatible with the one 
measured with the 2010 sample within statistical and systematic uncertainties, considering that the pp 
reference uncertainties are the same in the two measurements. In addition, the larger sample of central 
Pb-Pb collisions used in this analysis, compared to that used in the previous publication, enables the 
measurement of the D-meson Raa in a wider px range (the intervals 1 < px < 2 GeV/c and px > 
16 GeV/c were not accessible with the previous sample), with a substantial reduction (by a factor of 
about 2-3) of the statistical uncertainties. 


Figure |9] (left) also shows the average D-meson nuclear modification factor measured in minimum-bias 
p-Pb collisions at = 5.02 TeV 1133n . Since no significant modification of the D-meson production 
is observed in p-Pb collisions for px > 2 GeV/c, the strong suppression of the D-meson yields for 
Px > 3 GeV/c observed in central and semi-central Pb-Pb collisions cannot be explained in terms of 
cold nuclear matter effects and is predominantly due to final-state effects induced by the hot and dense 
medium created in the collisions. 


5.2 Comparison with results at lower collision energy 

In the right panel of Fig. |9l the average D-meson Raa for the 10% most central Pb-Pb collisions is 
compared to the D° nuclear modification factor measured by the STAR Collaboration for the 10% most 
central Au-Au collisions at = 200 GeV Bin . The D-meson Raa measured at the two energies 
are compatible within uncertainties for px > 2 GeV/c. It should be noted that the similar Raa of D 
mesons with high momentum, px > 5 GeV/c, i.e. in the range where the nuclear modification factor is 
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expected to be dominated by the effect of in-medium parton energy loss, does not necessarily imply a 
similar charm-quark energy loss at the two collision energies. Since the nuclear modification factor is 
also sensitive to the slope of the pj spectra in pp collisions, the combined effect of a denser medium and 
of the hai'der ^spectra at the LHC could result in similar values of Raa as at lower collision energies 
(see e.g. Ref. |^]). 


At low momentum (1 < px < 2 GeV/c), the Raa measured by STAR shows a maximum. This effect 
can be described by models including parton energy loss, collective radial flow and the contribution of 
the recombination mechanism to charm-quark hadronisation 0]. The ALICE results at higher do 
not show a maximum. However, the large uncertainties and the coarser binning at low px prevent a firm 
conclusion from being drawn. A different pattern could be explained by the different role of initial-state 
effects or of radial flow at the two collision energies. In the initial state, the modification of the parton 
distribution functions in a nuclear environment is predicted to lead to a stronger suppression of the heavy- 
quark production yields at low px with increasing 1641] . because of the smaller values of Bjorken-x 
probed. In addition, the momentum (kx) broadening effect, which gives rise to an enhancement of the 
R aa at intermediate px (Cronin peak), is known to be more pronounced at lower collision energies 11651 
661] . In the final state, in addition to energy loss, the collective expansion of the medium is also predicted 
to affect the momentum distribution of chai'med hadrons in heavy-ion collisions. Indeed, the interactions 
with the medium constituents are expected to transfer momentum to low-px charm quarks, which could 
take part in the collective radial flow of the medium. This effect could be enhanced by hadronisation 
via recombination, which is predicted in some models to contribute significantly to hadron formation 
at low and intermediate px OlSll . The momentum distributions of identified light-flavour hadrons at the 
LHC 11671. l68n indicate that the radial flow of the medium at LHC energies is about 10% higher than 
at RHIC 1691] . However, this stronger radial flow does not necessarily give rise to a more pronounced 
bump-like structure in the Raa at low px with increasing collision energy, because its effect can be 
counterbalanced by the different shape of the momentum spectra in pp collisions at different 

5.3 Comparison with pion and charged-hadron Raa 

As described in Section[TJ the colour-charge and quark-mass dependence of the energy loss can be tested 
with the comparison of D-meson and pion nuclear modification factors. In the left panel of Lig. [TOl 
the D-meson Raa (average of D*’, D+ and D*+) measured for the 10% most central Pb-Pb collisions is 
compared with the pion Raa in the interval 1 < px < 20 GeV/c and with the Raa of charged particles 
in 16 < Px < 40 GeV/c. The charged-particle Raa is shown in order to extend the comparison up to the 
higher px interval in which the D-meson yield was measured. The comparison of D mesons with charged 
hadrons at high-px is relevant because the Raa of different light-flavour hadron species are consistent 
with one another for px > 8 GeV/c 17211 . Moreover the contribution of pions dominates the charged- 
hadron yields at px of about 20 GeV/c with respect to other hadron species (about 65%) 17411 . A similar 
comparison is performed in the right panel of Pig.[T0]for the 30-50% centrality class. 

The Raa of D mesons and light-flavour hadrons are consistent for px > 6 GeV/c for both centrality 
classes. Lor px < 6 GeV/c, the Raa of D mesons tends to be slightly higher than that of pions. This can 
be also observed from the ratio of nuclear modification factors, presented in Pig. [T3j Considering that 
the systematic uncertainties of D-meson yields are mainly correlated with px, we observe 7?^^ > ^AA^t 
low Px with a significance of about 1 a in four px intervals, in the most central events. In the 30-50% 
centrality class, the significance of the effect is smaller than in central collisions. 

A direct interpretation of a possible difference between the D-meson and pion Raa at low px is not 
straightforward. In the presence of a colour-charge and quark-mass dependent energy loss, the harder 
Px distribution and the harder fragmentation function of charm quarks compared to those of light quarks 
and gluons could lead to similar values of D-meson and pion Raa, as discussed in Ref. 11311 . In addition, 
it should be considered that the pion yield could have a substantial contribution from soft production pro- 
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Figure 9: Left; prompt D-meson /?aa (average of D*', D+ and D*+) as a function of pj in Pb-Pb collisions 
at y'sNN = 2.76 TeV in the 0-10% and 30-50% centrality classes. Prompt D-meson nuclear modification factor 
(average of D**, D+ and D*+) as a function of pi in p-Pb collisions at 


= 5.02 TeV 113311 . Right: prompt 
D-meson Raa (average of D**, D+ and D*+) as a function of pi in the 10% most central Pb-Pb collisions at 
= 2.76 TeV compared to D** Raa measured by the STAR Collaboration in Au-Au collisions at RHIC at 
i/iNN = 200 GeV im. A zoomed-in plot of the interval 0 < px < 8 GeV/c is shown in the inset. Statistical (bars), 
systematic (empty boxes), and normalisation (shaded boxes at Raa =1) uncertainties are shown. Horizontal bars 
represent bin widths. Symbols are placed at the centre of the bin. 




Figure 10: Prompt D-meson Raa (average of D*', D+ and D*+) as a function of pi compared to the nuclear 
modification factors of pions 172 ] and charged particles iF^ in the 0-10% (left) and 30-50% (right) centrality 
classes. Statistical (bars), systematic (empty boxes), and normalisation (shaded box at Raa = 1) uncertainties are 
shown. Horizontal bars represent bin widths. Symbols are placed at the centre of the bin. 


cesses up to transverse momenta of about 2-3 GeV/c due to the strong radial flow at LHC energies. This 
soft contribution, which is not present in the D-meson yield, does not scale with the number of binary 
nucleon-nucleon collisions. Finally, the effects of radial flow and hadronisation via recombination, as 
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Figure 11: Average of prompt D**, D+ and D*+ f?AA in the centrality classes 0-10% (a and b) and 30-50% (c 
and d) compared with model calculations: Djordjevic (75|, CUJET3.0 ||^, 77], WHDG Vitev if^ (a and 

c), TAMU elastic S, Cao, Qin and Bass S, MC@sHQ+EPOS, Coll+Rad(LPM) S, POWLANG [ill HI, 
BAMPS PHSD (b and d). Some of the model calculations are shown by two lines to represent their 

uncertainties. 


well as initial-state effects, could affect D-meson and pion (light-flavour particle) yields differently at a 
given pj, thus introducing an additional complication in interpreting the magnitude of the Raa in terms 
of different in-medium parton energy loss of charm quarks, light quarks and gluons. 


5.4 Comparison with models 


Figure [TT] shows the comparison of the average D-meson Raa for the two centrality classes 0-10% (a 
and b) and 30-50% (c and d) with most of the available model calculations. The model calculations are 
described and compared in a recent review lllSn . A concise summary is given in the following paragraphs. 


The interaction of heavy quarks with the medium constituents is computed considering radiative and 
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collisional processes in the calculations indicated as Djordjevic [TJ], WHDG i20l422ll . CUJET3.0 11761 
Wh . MC@ sHQ+EPOS BAMPS | 53 - 3 ]> and Cao^Q^, Bass |79]. Only collisional interactions 
are considered in the model calculations POWLANG |^, 82], TAMU elastic |70] and PHSD |83]. In 
BAMPS, two different options are considered: including only collisional energy loss but introducing a 
scaling factor to match RHIC high-pT data (where radiative energy loss is expected to be dominant) 
or including both collisional and radiative energy loss. Also for the Vitev model iVSH two different 
options are considered: including only radiative energy loss {Vitev rad) or also considering the in-medium 
dissociation of heavy-flavour hadrons {Vitev rad-i-dissoc). 

The medium is described using an underlying hydrodynamical model in CUJET3.0, Cao, Qin, Bass, 
MC@sHQ+EPOS, BAMPS, POWLANG, TAMU elastic and PHSD, while Djordjevic, WHDG and Vitev 
use a Glauber model nuclear overlap without radial expansion. 

The initial heavy-quark pj distributions are based on next-to-leading order (NLO) or FONLL pertur¬ 
bative QCD calculations in all model calculations, except for Cao, Qin, Bass, which uses the PYTHIA 
event generator The EPS09 NEO parameterisation ll^ of the nuclear parton distribution functions 
is included by POWLANG, MC@sHQ+EPOS, TAMU elastic, PHSD and Cao, Qin, Bass. 

All model calculations use in-vacuum fragmentation of heavy quarks for the high-momentum region. 
At low momentum this is supplemented by hadronisation via recombination in the MC@sHQ+EPOS, 
POWLANG^ Cao, Qin, Bass, TAMU elastic and PHSD models. The two last models also include 
scattering of D mesons in the hadronic phase of the medium. Also for the Cao, Qin, Bass model, the 
hadronic-rescattering effects have been studied in a recent publication 18411 and no large differences in 
the /?aa are observed, when these processes are considered. 


Several model calculations provide a good description of the measured Raa for both centrality classes. 
The MC@sHQ+EPOS model has recently improved the description of the /?aa in the pj interval 2- 
8 GeV/c including the EPS09 shadowing parameterisation in addition to in-medium energy loss, the 
TAMU elastic model overestimates the /?aa in central collisions in the pi interval 6-30 GeV/c and the 
POWLANG model underestimates it in the interval 5-36 (8-16) GeV/c in the 0-10% (30-50%) centrality 
class. Interestingly, these model calculations provide a fair description of the D-meson V 2 measured at 
EHC i35n and of the D-meson /?aa measured at RHIC OlH . On the other hand, the model calculations 
that do not include a hydrodynamical medium expansion and hadronisation via recombination, namely 
Djordjevic, Vitev, WHDG —and as a consequence do not describe the features observed for the V 2 at the 
EHC and the Raa at RHIC in the momentum region up to about 3-5 GeV/ c — provide a good description 
of the /?aa in the full “high pj interval”, above 5 GeV/c. The Vitev model shows a better agreement when 
including the D-meson in-medium dissociation mechanism. The BAMPS model with collisional energy 
loss describes the data better for the low-pT interval, as is the case for the D-meson V 2 11351]. The inclusion 
of radiative energy loss improves the agreement at high pj. The Cao, Qin, Bass model describes the Raa 
in both centrality classes, but underestimates the D-meson V 2 1350. The PHSD model describes the Raa 
in both centrality classes. 


Figure[T2]shows the TAMU elastic and MC@sHQ+EPOS calculations of the nuclear modification factor, 
for the 10% most central events, with and without including the EPS09 shadowing parameterisation. 
For both models the inclusion of shadowing reduces the /?aa by up to about 30-40% in the interval 
Pt: <5 GeV/c, resulting in a better description of the data. 


Four of the model calculations also provide the nuclear modification factor of pions and charged particles 
{Djordjevic, CUJET3.0, WHDG and Vitev). All these calculations include radiative and collisional energy 
los^ The left panel of Fig. [13] shows the comparison with the measured charged-pion 7?^^ (pj < 

^Note that recombination was not included in the version of the POWLANG model used for the comparison with the D- 
meson V2 measurement in j^. 

^The in-medium formation and dissociation process, included by Vitev for D mesons, is not relevant for pions, which have 
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Figure 12: Average of prompt D**, D+ and D*+ f?AA in the centrality classes 0-10% compared with TAMU elastic 
and MC@sHQ+EPOS models calculations with and without including EPS09 shadowing parameterisations 16411 . 




Figure 13: Left: /?aa of charged pions {pi < 16 GeV/c) iF^ and of charged particles (px > 16 GeV/c) iF^ 
compared with model calculations that compute also the D-meson /?aa- Right: ratio of the Raa of prompt D 
mesons (average of D**, D+ and D*+ as shown in Fig. fTTIi and the Raa of charged pions (for px < 20 GeV/c) or 
charged particles (for px > 20 GeV/c), compared with the same model calculations shown in the left panel. 


16 GeV/c) 11721] and charged-particle 7 ?aa (Pt >16 GeV/c) 1173h . The model calculations provide a 
reasonable description of the measurements, with WHDG generally showing smaller /?aa values than 
seen in data, although consistent within experimental and theoretical uncertainties. 


The right panel of Fig. [13] shows the ^aa/^aa FPt < 16 GeV/c) and /?aa/^aa (Pt > 16 GeV/c) ratios 
for data and for these four model calculations. In the case of data, the uncertainties of D-meson and 
charged-pion (or charged-particle) measurements were propagated as uncorrelated uncertainties, except 


a much larger formation time. 
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for the uncertainty on (Taa), which cancels in the ratic0. In the case of model calculations, the theoretical 
uncertainty, when provided, was propagated assuming full correlation between D mesons and pions 
(charged particles), since it accounts for a variation of the medium density (or temperature). Only the 
Djordjevic and CUJET3.0 models, which use radiative and collisional energy loss, can describe the two 
Raa results and their ratio over the full pj interval in which they provide the calculations (pj > 5 and 
8 GeV/c, respectively). The Vitev model can describe the data at the lowest pj (2-6 GeV/c) only if the 
dissociation mechanism is included, suggesting that the effect is relevant in this model. However, the 
model overestimates the data in the interval 6-12 GeV/c. 


6 Conclusions 


We have presented the measurements of the production of prompt D**, D+ and D*+ mesons at central 
rapidity in Pb-Pb collisions at a centre-of-mass energy per nucleon pair = 2.76 TeV, as well as 
their nuclear modification factor Raa- The measurements cover the interval 1 < pt < 36 GeV/c for the 
0-10% centrality class and 1 < pj < 16 GeV/c for the 30-50% centrality class. 


The nuclear modification factor shows a maximum reduction of the yields with respect to binary scaling 
by a factor 5-6, for transverse momenta of about 10 GeV/c for the 10% most central Pb-Pb collisions. 
A suppression of a factor about 2-3 persists in the highest px interval covered by the measurements 
(24-36 GeV/c). At low px (1-3 GeV/c), the Raa has large uncertainties, that span the range from 0.35 
(factor of three suppression) to 1 (no suppression). In all px intervals above 5 GeV/c, the Raa for the 
30-50% centrality class is about twice that for the 0-10% centrality class. The suppression observed for 
Px > 3 GeV/c is interpreted to be due to interactions of the charm quarks within the high-energy density 
medium formed in the final-state of Pb-Pb collisions. This is demonstrated by the nuclear modification 
factor measurements in p-Pb collisions at = 5.02 TeV, which indicate that D-meson production is 
consistent with binary collision scaling ll33n . 


The D-meson Raa was compared with that of charged pions in the interval 1 < px < 16 GeV/c, also 
in terms of the ratio /?aa/^aa’ charged particles up to px = 36 GeV/c (7?^^//?^^). In 

the interval 1 < px < 6 GeV/c, the Raa values of D mesons are higher than those of pions, although 
consistent within uncertainties. For the 10% most central collisions, the ratio R^a/^aa larger than 
unity by about 1 a of the total uncertainties, which are to some extent correlated among px intervals. 
For Px > 8 GeV/c, the Raa values are compatible with those of pions and charged particles up to 
Px = 36 GeV/c. 


Several models provide a good description of the Raa for both centrality classes. Interestingly, the 
models that show larger deviation from the data, especially in the high-px region, are among those that 
provide a good description of the D-meson V 2 measured at the FHC and of the D-meson Raa measured 
at RHIC, in the low-px region. On the other hand, the models that do not include a hydrodynamical 
medium expansion and recombination, and as a consequence do not describe V 2 in the momentum region 
up to about 3-5 GeV/c, provide a good description of the Raa at the FHC in the full high-px interval, 
above 5 GeV/c. 

Only two out of the four models that compute R^a/^aa describe this measurement over the full px 
interval for which they provide the calculations. In these models, the nuclear modification factors of D 
mesons and pions turn out to be very similar as a consequence of a compensation among the larger energy 
loss of gluons with respect to that of charm quarks (mainly due to the larger colour coupling factor), the 
different amount of gluon and light quark yields on the pion Raa and the harder px distribution and 
fragmentation of charm quarks with respect to those of gluons and of light quarks. 


^The uncertainty on the normalisation (integrated luminosity) of the pp reference cross sections for D mesons and pions 
(charged particles) does not cancel in the ratio, because the two cross sections were measured in two data samples at different 
centre-of-mass energies. 
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